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Abstract 


A time-iterative full Navier-Stokes code, PARC, is used to analyze the 
flowfield of a two-dimensional ejector nozzle system. A parametric study was 
performed for two controlling parameters, duct to nozzle area ratio and nozzle 
pressure ratio. Results show that there is an optimum area ratio for the effi- 
cient pumping of secondary flow. At high area ratios, a freestream flow passes 
directly through the mixing duct without giving adequate pumping. At low area 
ratios, the jet boundary blocks the incoming flow. The nozzle pressure ratio 
variation shows that the pumping rate increases as the pressure ratio increases, 
provided there is no interaction between the shroud wall and the shock cell 
structure . 


Introduction 


A steady increase in market needs for trans-Pacific aircraft operation 
renewed national interest in developing a high-speed civil transport (HSCT) . 

Two of the main environmental concerns are airport community noise and engine 
emission (ozone depletion). This study is related to the community noise. 

The acoustic challenge to designing HSCT engines is to reduce noise to FAR 36 
stage III noise levels while maintaining the economic viability of the air- 
craft. To meet these requirements and still provide a high level of perform- 
ance, the development of a low-noise exhaust nozzle design is essential. To 
date, various noise suppression concepts have been investigated but none have 
reduced the noise to the FAR 36 stage III requirement in an efficient manner. 
Progress and development of these technologies are well reviewed in Ref. 1. 
However, one of the more current, promising technologies is the mixer-ejector 
concept developed by W. Presz et al.2 - 4 The mixer-ejector concept introduces 
an array of large-scale, low-intensity streamwise vortices into the downstream 
mixing duct. These vortices enhance mixing through a convective stirring proc- 
ess. This process results in increased pumping performance with more com- 
pletely mixed flows exiting the ejector shroud, which consequently reduces jet 
exhaust noise. Various experimental studies are presently being performed by 
the NASA Lewis Research Center and others to verify this concept and to make it 
practical for application to HSCT. 

Numerous degrees of freedom as well as complexity in the experiment, how- 
ever, make it difficult to test every possible set of combinations. An analyt- 
ical or computational approach is necessary to aid the understanding of 
experimental results and to provide some guidelines for the ongoing experiment. 
In the present study as a preliminary step to analyzing the flowfield inside 
the mixer-ejector nozzle, the simplified two-dimensional ejector tested at NASA 
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Lewis is considered first. A schematic of this two-dimensional ejector nozzle 
is shown in Fig. 1. The ejector is composed of three major components - pri- 
mary nozzle, secondary inlet, and mixing duct. The secondary flow is drawn 
from the freestream. The inlet has an elliptical bellmouth to direct the sec- 
ondary airflow smoothly into the mixing duct. This mixing duct has a diffuser 
section which lowers the mixing duct pressure to enhance the ejector pumping 
rate. The diffuser wall angle for the present case is 4°. The mixing duct 
area can be varied by moving the shroud vertically. 

The flow feature in the ejector nozzle system is complicated by the strong 
inviscid/viscous interaction, freestream entrainment, shock cell structure, and 
the coexistence of supersonic and subsonic regions. Until now, most of the 
analytical studies for ejector nozzle systems were based on Von Karman's control 
volume approach. ^>6 It is necessary, however, to use the complete Navier-Stokes 
equations to capture the previously mentioned complex phenomena properly. 
Recently, Deese and Agarwal 7 presented numerical solutions of Navier-Stokes 
equations for ejector flowfields using a Runge-Kutta time stepping method; 
their results were for the subsonic flow regime. In the present study, a time- 
iterative full Navier-Stokes equation solver (PARC code) is used to investigate 
the ejector flowfields for mainly underexpanded flows. 

Features of the ejector flow depend on several parameters - e.g., geometry 
of the mixing duct, duct to nozzle area ratio, duct length to width ratio, pri- 
mary nozzle location with respect to the inlet of mixing duct, primary jet 
pressure, and Reynolds number, etc. In this study, the effects of the two con- 
trolling parameters duct to nozzle area ratio (AR) and nozzle pressure ratio 
(NPR) on the ejector flowfield were investigated. Computational results are 
compared with experimental results to validate the code. 


Governing Equations 

The flow through the ejector nozzle can be described by the Navier-Stokes 
equations. The nondimensional form of the Navier-Stokes equations can be writ- 
ten in the strong conservation-law form and in the generalized coordinate sys- 
tem as follows: 
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where the primary dependent vector Q is given by 
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The vectors E and F represent inviscid flux vectors: 
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In these expressions, x and y are the Cartesian coordinates, and £ 
and n are the transformed coordinates. The quantity J is the Jacobian of 
the transformation defined as 3(£ ,n)/3(x,y) . The variables p and p 
denote density and static pressure, u and v are velocity components in 
Cartesian coordinates, and e is the total energy per volume defined to be 
p[i + (u 2 + v 2 )/2] , where i is the specific internal energy. The contravar- 
iant velocities U and V are defined as 


U = C x u + £yV V = Hx u + HyV (4) 

On the right side of Eq. (1) the vectors and V n represent viscous 

flux vectors composed of linear combinations of V x and Vy which are the 
Cartesian components of the viscous flux vectors: 

V=-(iV+iv') V -IfriV + iv') (5) 
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The Cartesian viscous vectors V x and Vy are 



where the components of the shear stress tensor are given by 
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Here p is the viscosity, k is the thermal conductivity, and X is the sec- 
ond coefficient of viscosity. The Re and Pr are the Reynolds number and 
Prandtl number, respectively. For the present study, the property p is eval- 
uated by the Sutherland formula as 
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where T s is the Sutherland temperature constant. The Navier-Stokes equations 
are augmented by the equation of state taken here as that of a perfect gas. 
Thus, the total energy is defined by 


e = 



(9) 


where y is the ratio of specific heats. 


In the turbulence calculations, the algebraic turbulence model based on 
the Thomas formulation of the Baldwin and Lomax model®’9 is employed for eddy 
viscosity. 


Numerical Procedure 


The numerical solution of the Navier-Stokes equations is obtained by the 
PARC computer code, which is a variant of the ARC code developed by Pulliam 
et al.10 The detailed development and some of the recent work related to the 
PARC code can be found in Ref. 11. In general, the program uses an Euler 
implicit discretization in time along with central differencing in space for 
both the inviscid and viscous terms. An efficient solution of the resulting 
matrix can be obtained by using the approximate factorization based on the 
Douglas-Gunn procedure 13 with diagonal izat ion of the inviscid terms. In addi- 
tion, Jameson's type of artificial dissipation 13 is added and treated impli- 
citly for monotonicity and stability. This procedure results in the scalar 
pentadiagonal matrix which can be inverted in an efficient manner. All equa- 
tions are solved simultaneously at each time step until convergence to a steady 
state is reached. 


Boundary Conditions 

The flow domain of interest is shown in Fig. 1. In the present study, we 
compute only half of this domain since the flowfield is symmetric. Boundary 
conditions are specified at the upstream and downstream ends, at the far field, 
and on all walls. At the upstream boundary stagnation pressure, the stagnation 
temperature and flow angle (v/u) of the incoming stream are specified. These 
three boundary conditions are supplemented by extrapolating the left running 
Riemann invariant. The flowfield at the downstream end is mixed subsonic/ 
supersonic. For the supersonic flow region, since no physical boundary condi- 
tion is needed, all flow variables are extrapolated. For the subsonic flow 
region, back pressure is specified for a given freestream Mach number. The no 
slip condition and thermally adiabatic conditions are imposed on all walls. 

This condition is again augmented by the normal momentum equation. At the noz- 
zle exit and at the trailing edge of the shroud, a Kutta condition is imposed 
to ensure the continuity of pressure on both sides of the wall: 
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The far-field boundary condition for the present case uses normal velocity 
v = 0, and the remaining flow variables are extrapolated. A symmetry condition 
is imposed on the plane of symmetry of the computational domain. All the 
boundary conditions mentioned previously are imposed in an explicit manner. 
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Results and Discussion 


To understand the physical flow phenomena in the ejector nozzle system, a 
parametric study was done for two controlling parameters, mixing duct to nozzle 
area ratio (AR) and nozzle pressure ratios (NPR). Herein, AR is defined as the 
ratio of the mixing duct inlet area (As) to the primary nozzle throat area (Ap) 
(AR = As/Ap, refer to Fig. 1). All calculations are performed on a 171 x 91 
grid. The typical grid used for the present calculations is shown in Fig. 2. 
Grid is clustered near all walls (the first grid line is located at a y + of 
approximately 10.0) and at the nozzle exit to resolve boundary layer and rapid 
expansion of the flow. The far-field boundary is located about five widths of 
the mixing duct from the shroud so that wall effects may be negligible. In 
all calculations, a freestream Mach number of 0.2 is used to simulate takeoff 
conditions and the flow angle at the upstream end is v/u = 0. 

Figure 3 shows results for AR's ranging from 1 to 6. Varying the area 
ratio was accomplished by moving the shroud wall vertically relative to the 
centerline. Figure 3 shows Mach numbers and pressure contours for the differ- 
ent area ratios at a nozzle stagnation pressure of 3.5 atm and an ambient 
stagnation pressure of 1 atm, which corresponds to an NPR of 3.5. The nozzle 
stagnation temperature is 505 K, and the ambient stagnation temperature is 
300 K. In all calculations, the Reynolds, Prandtl, and turbulent Prandtl num- 
bers used are 1.3xl0 7 , 0.72, and 0.9, respectively. The Reynolds number is 
based on the nozzle width and nozzle exit velocity which would be obtained for 
an i sen tropic expansion. 

For all area ratios, the Mach number and pressure contours inside the pri- 
mary nozzle remain unchanged due to the choking at the exit. For the present 
calculations, the nozzle pressure ratio exceeds the critical pressure ratio 
(1.9); thus, the flow beyond the nozzle exit becomes supersonic and results in 
the formation of a series of shock cells. In the present case, the jet bound- 
ary is that of a constant pressure and the waves reflected on this boundary 
become waves of the opposite senses (compression waves are reflected as expan- 
sion waves). This mechanism forms a shock cell structure in the underexpanded 
flow. At AR = 6, the primary flow expands rapidly at the nozzle exit and 
becomes a jet with a series of shock cell structures. At the secondary flow 
passage, however, the Mach number is nearly uniform inside the mixing duct. 

This suggests that freestream flow is passing directly through the mixing duct; 
hence, mixing between the primary and secondary flows is quite small. The 
pressure contour also shows nearly uniform pressure through the secondary flow 
passage. This again indicates uniform flow is passing through without giving 
efficient pumping. 

The results at AR = 2.5 are similar to those at AR = 6 . The Mach number 
contours in the mixing duct still show the uniform flow region, but the incom- 
ing flow is more accelerated due to the effective reduction in the secondary 
flow area between the shroud wall and the jet boundary. Consequently, the rate 
of mixing increases. The pressure contour also shows an enlarged low pressure 
region in the mixing duct inlet region, suggesting more effective pumping. 

At AR = 1.5, the Mach number changes continuously inside the mixing duct 
with much faster incoming flow at the mixing duct inlet. Also, the pressure 
contours show that the mixing duct inlet region is dominated by low pressure, 
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which leads to more efficient pumping. However, there is a small separation 
bubble at the trailing edge region which may deteriorate the thrust perform- 
ance. This separation occurs because of a very strong adverse pressure gradi- 
ent (shown later in Fig. 5). 

The Mach number contour at AR = 1 shows that the jet boundary nearly 
impinges on the shroud wall. This behavior reduces the effective secondary 
flow area much less than previous cases and the Mach number increases abruptly. 
However, Fig. 5 shows that the Mach number at the mixing duct inlet is smaller 
than that when AR = 1.5. This suggests that the pumping effect is also 
smaller for AR = 1 than it is for AR = 1.5. The reason may be that the jet 
boundary, which is hitting the shroud wall, blocks the incoming flow. Also, 
the shock cell structure is significantly more distorted than those at larger 
area ratios. Unlike the previous cases, the jet boundary is no longer at a 
constant pressure, and the pressure varies continuously inside the mixing duct 
as shown in Fig. 3. This causes a significant change in the shock cell struc- 
ture. Pressure contours also show that pressure at the mixing duct inlet is 
higher at AR » 1 than it is at AR = 1.5, which again suggests less effective 
pumping at AR = 1 . 

Some representative streamwise velocity profiles for these cases are shown 
in Fig. 4. For the four area ratio cases, velocity profiles are presented at 
the three different axial locations noted in Fig. 1. The profile at station 1 
is at the exit of the convergent nozzle, and the profiles of stations 2 and 3 
are inside the mixing duct region. The streamwise velocity profiles, again, 
demonstrate the relatively slow rate of mixing between the primary and second- 
ary flows at high area ratios and the increasingly rapid development of flow at 
the low area ratios. When the velocity magnitudes at the mixing duct inlet 
section are compared, the average velocity at AR = 1.5 is higher than those 
velocities at other area ratios. These findings indicate that AR = 1.5 is at 
the optimum shroud position for the pumping capability at this NPR. 

The calculated shroud surface pressure distribution is shown in Fig, 5. 
Cases for four different area ratios have been presented. At AR = 6 and 2.5, 
the pressure decreases with downstream distance to reach the minimum near the 
inlet section of mixing duct; this causes the secondary flow to accelerate in 
the inlet region. The surface pressure then increases as the flow proceeds 
downstream until it reaches ambient pressure at the ejector outlet. At 
AR = 1.5, the pressure distribution along the shroud wall is quite different 
from that for other area ratios. A large pressure dip is found at about 
30-percent chord length of the shroud after a smooth pressure decrease in the 
inlet region. At these conditions, the interaction between the shock cell 
structure of the jet and the shroud wall predominates, which further reduces 
the secondary flow passage. The pressure dip is attributed to the previous 
interaction phenomena. At AR = 1 , the pressure dip is remarkable. The jet 
is nearly hitting the shroud wall (refer to Fig. 3). This very strong interac- 
tion between the jet and the shroud wall greatly decreases the secondary flow 
passage. Flow acceleration through this passage causes the pressure to drop 
very quickly (near to the sonic condition). Also, the pressure in the mixing 
duct inlet is higher than it is for the AR - 2.5 and 1.5 cases, suggesting a 
poor pumping characteristic. The pressure distribution along the airfoil 
shroud wall is compared with the experimental data obtained in the NASA Lewis 
9 by 15 foot low-speed wind tunnel at corresponding area and nozzle pressure 
ratios (Fig. 5). The agreement is quite good at AR = 2.5 over most of the 
shroud wall. 
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A similar set of calculations is shown in Fig. 6 for the case where the 
nozzle pressure ratio is varied at a constant area ratio. Four nozzle pressure 
ratio cases have been computed for NPR = 1.5, 2.5, 3.5 and 4.5. Corresponding 
nozzle stagnation pressures are 1.5, 2.5, 3.5 and 4.5 atm, while the ambient 
stagnation pressure is kept as 1 atm. In all calculations, AR = 2.5 is used. 

Comparing Mach number contours at various pressure ratios shows that as 
the nozzle pressure increases the shock cell structure becomes larger. At 
NPR =1.5, the shock cell structure does not appear because this pressure 
ratio is below the critical pressure and the jet is entirely subsonic. At 
NPR = 2.5, the shock cell begins to show up in the vicinity of the nozzle exit; 
at higher pressure ratios, 3.5 and 4.5, the shock cell becomes larger and prop- 
agates further into the downstream region. As mentioned before, this large 
shock cell structure causes a reduction in the secondary flow area between the 
shroud wall and the jet boundary and, therefore, the incoming flow into the 
mixing duct accelerates as the pressure ratio increases. This suggests more 
efficient pumping at higher pressure ratios. Also, the jet spreading rate 
increases as the pressure ratio increases. 

Pressure contours are also shown in Fig. 6. The development of a shock 
cell structure can be seen more clearly in these pressure contours. At 
NPR =1.5, pressure is nearly uniform inside the mixing duct region without any 
shock cell. But as the pressure ratio increases, the shock cell becomes larger 
and stronger. In general, as the pressure ratio increases, the pressure at the 
mixing duct inlet region decreases. This condition leads to more efficient 
pumping at higher pressure ratios. 

The corresponding pressure distribution on the shroud surface is shown in 
Fig. 7. In general, the pressure at the mixing duct inlet region decreases as 
the nozzle pressure ratio increases. As mentioned before, this indicates more 
pumping at higher pressure ratios. At NPR = 4.5, however, we see a similar 
behavior to what was seen at AR = 1.5 in Fig. 5. The pressure is at a minimum 
at around 30 percent of shroud chord length. This occurs because at a high 
pressure ratio the strong shock cell structure decreases the secondary flow 
passage and hence the flow accelerates. 

The pumping characteristic is an important design parameter of the ejector 
nozzle for jet noise reduction as well as for ejector performance. Figure 8 
shows the pumping characteristic of the ejector nozzle system for various area 
ratios and nozzle pressure ratios. Computations were made for AR's between 1 
and 12 and for NPR's of 2.5, 3.5 and 4.5. The pumping characteristic is 
expressed in terms of entrained mass flow per unit area (Ws/As). This parame- 
ter is a better indicator of pumping than the secondary to primary mass flow 
ratio (Ws/Wp) when freestream is present. The pumping characteristic is quite 
low and approaches a certain asymptotic value at high area ratios. This is 
because at high area ratios the pumping is negligible and most of the mass 
flow in the secondary passage is due to freestream flow. But as the area 
ratio decreases, pumping increases and reaches a peak value, depending on noz- 
zle pressure ratio. (At NPR =3.5 the peak pumping rate occurs at about 
AR = 1.5; at NPR =4.5 the peak pumping rate occurs at about AR = 2.5; and 
at NPR =2.5, the pumping rate increases until AR = 1 is reached.) A poor 
pumping characteristic is seen for very low area ratios at NPR =3.5 and 4.5 
because of the blocking effect caused by the shroud wall and the shock cell 


7 



structure interaction. It is also anticipated that there is a certain area 
ratio for NPR =2.5 where pumping will be reduced. 


Summary 

A computational study of the flowfield in a two-dimensional ejector nozzle 
has been described. The flow feature inside the ejector nozzle system is com- 
plicated by strong inviscid/viscous interaction, freestream entrainment, shock 
cell structure, and the presence of both supersonic and subsonic flow regions. 

A time-iterative full Navier-Stokes code, PARC, is used to capture this complex 
phenomena. A parametric study was made for two controlling parameters, duct to 
nozzle area ratio and nozzle pressure ratio. Results show that there is an 
optimum area ratio for the efficient pumping of secondary flow. At high area 
ratios, the freestream flow passes directly through the mixing duct without 
providing adequate pumping. At low area ratios, the jet boundary blocks incom- 
ing flow. The nozzle pressure ratio variation shows that the pumping rate 
increases as the pressure ratio increases - provided there is no interaction 
between the shroud wall and the shock cell structure (i.e., for a certain area 
ratio). The comparison with experimental data shows quite good agreement. 
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Fig. 2. Computational grid for ejector nozzle (171 x 91). 
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Mach number 


Pressure contours 




(c) AR = 1.5. 





(d)AR- 1.0. 

Fig. 3. Mach number and pressure contours for area ratios (AR) of 6, 2.5, 1.5, and 1 at a nozzle pressure ratio of 3.5. 
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Mach number 


Pressure contours 



(b) NPR - 2.5. 



(c) NPR - 3.5. 




(d) NPR - 4.5. 


Fig. 6. Mach number and pressure contours for nozzle pressure ratios (NPR) of 1.5, 2.5, 3.5, and 4.5 at an area ratio of 2.5. 
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Fig. 7. Pressure distribution on the shroud surface. 


170 

□ 


160 

< D 

1 


150 

- * O 

o A 


140 

A 


130 

□ 

o 

J 

A 

120 


□ 

110 

O 


100 

90 

° 1 



90 I 1 L 

0 2 4 


Nozzle 

pressure 

ratio, 

O 4.5 

A 3.5 

□ 2.5 


0 

A 

□ 

8 

□ 

1 J 1 

6 8 10 

Area ratio 


6 

J 

12 


Fig. 8. Pumping characteristics of ejector nozzle at 
various area and nozzle pressure ratios. 
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